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We propose a feasible scheme to achieve holonomic quantum computation in a decoherence-free 
subspace (DFS) with trapped ions. By the application of appropriate bichromatic laser fields on the 
designated ions, we are able to construct two noncommutable single-qubit gates and one controUed- 
phase gate using the holonomic scenario in the encoded DFS. 
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Realization of practical quantum computers presents a 
tremendous challenge to scientists. One of the most dif- 
ficulties is caused by the random errors during gate op- 
erations, which mainly roots in decoherence and stochas- 
tic errors in the control process. To suppress the latter, 
the holonomic quantum computation(HQC) P| as well as 
a more generalgcomctric quantum computation(GQC) 
scenario 0i S 3 were proposed, which arc respectively 
based on the adiabatic non-Abelian holonomy and the 
nonadiabatic geometric phase. In the case of HQC, to 
realize the wanted gate operation, the Hamiltonian is 
driven adiabatically to undergo a designated cycle in the 
controllable parameter space, in which degenerate dark- 
states (with zero eigenvalues) may be defined. Thus for 
an initial state, only the adiabatic geometric phases can 
be accumulated in the designated operation, and the cor- 
responding holonomic quantum gates may be achieved. 
Because of its global geometric feature, it is believed that 
the HQC scheme is rather robust against the stochastic 
errors occurring in the adiabatic cyclic evolution^. On 
the other hand, decoherence is another main obstacle for 
quantum computation. To keep a quantum system away 
from the main decoherence source, one should attempt 
to isolate the system from the environment that causes 
mainly the decoherence effect. To deal with this mat- 
ter, several scenarios were put forward, including quan- 
tum error correcting codes [6j and quantum error avoid- 
ing codes. In particular, when the system-environment 
interaction has certain symmetry, a promising scheme 
based on the decoherence-free subspaces(DFS) was pro- 
posed 0, H, 0- Depending on the type of interaction 
symmetry, one sort of DFS is usually immune to certain 
system-bath disturbances. 

Very recently, combining the HQC with a four-qubit- 



encoding DFS and based on a kind of model Hamiltonian, 
Wu et al proposed an interesting strategy to realize uni- 
versal quantum computation Such a strategy was 
indicated to afford the ability against the collective de- 
coherence induced by environment and to be rather ro- 
bust to certain stochastic errors during operations. Mo- 
tivated by the above idea and considering that the four- 
qubit DFS scenario may hardly be scalable ^3 j we here 
propose a feasible scheme based on a simple two-qubit- 
encoding DFS and the adiabatic non-Abelian holonomy, 
which is distinctly different from the four-qubit-cncoding 
scheme. In particular, by the application of appropri- 
ate laser fields, wc derive an effective laser-ion Hamilto- 
nian in a certain kind of trapped ion systems and illus- 
trate how to physically implement universal honolomic 
quantum computation in the present two-qubit-encoding 
DFS. 
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FIG. 1: A schematic illustration of the encoded sub-systems 
(logical qubits) in trapped ions. Each sub-system consists 
of two ions enclosed by the dashed circle and the effective 
interaction between any two logical qubits is denoted by the 
dotted line. 

To construct the hololomic quantum gates using the 
decoherence-free subspaces scheme in a linear Paul trap, 
we need first to obtain the wanted effective Hamiltonian 
for a sub-system consisting of two ions (e.g., ions ji and 
j2 in Fig.l) with the appropriate laser fields on. Let us 
assume that each ion has the four-level energy structure, 
as shown in Fig. 2 [TH . The Hamiltonian may be written 
as [ll(?i= 1) 
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FIG. 2: Schematic level structure for single ions. |e), |0), and 
|1) as well as \a) denote the exited state, the lowest ground 
state, and the other two stable (metstable) states (with the 
last one as an ancillary state), respectively. [0) and |1) are 
the basic units used to encode the logical qubits of DPS. 
Sic (a = eO, el,ea) are the corresponding Rabi-frequencies 
with the resonant laser fields being applied on the ions. 
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Ho = i'{a+a+ 1/2) + ^Ea\a)i{a\ 



(2) 



H^nt = ^f^L(e'f"'^"+"^^""'*'|e).(a| + H.c), (3) 

where a € {0,1, a} denote the inner states of an ion, v 
is the frequency of the ccnter-of-mass vibrational mode 
of ions chain, a"*" and a are the creation and annihilation 
operators of the collective motion, Ea is the energy of the 
I Of) state, rigQ, is the resonant Rabi frequency between the 
states a and e of the ith ion in the laser field with the 
frequency Ui. The positions of the ions Xi are replaced 
by ladder operators 'r]i{a+a~^) |l2l |. where ry^ is the Lamb- 
Dicke parameter denoting the ratio of the ion oscillation 
and the wave length of the exciting radiation. Since we 
tune the lasers close to the the center-of-mass vibrational 
mode where all ions vibrate in the same way, the coupling 
to the vibration is uniform for all ions, i.e., rji = 77 for all 
i. 

We now choose the so-called bichromatic laser field ^3 
to couple the two ions in the same sub-system (see 
Fig. 2). For the jith ion, we let it be illuminated by 
two laser pulses, which are resonant to the transition be- 
tween |1) and |e), with frequency ujj-^ = tUei + — <5io), 
ujj_^ ~ ujei — {ly ~ Sio), where Wei denotes the energy dif- 
ference between the state |e) and |1) and 610 represents 
the additional detuning. In the Lamb-Dicke limit that 
ri'^{n+l) <C 1, where n is the vibrational numbers, we can 
approximately expand the exponential gene e*(''(°+° to 
1 -|- iri{a -I- a"*"). Under the rotating- wave approximation 
and in the interaction picture, the laser-ion interaction 
of the ji th takes the form 



+ a+e^'^-']\e),AM+H.c.. (4) 

At the same time, we also apply two laser pulses on the 
j2^/i ion with a;,2 = Ueo + [i^ - Sio), uj. = uJeo - - Sio), 



where uieo denotes the energy difference between the state 
|e) and |0). Similarly we can have 
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H.c. 



(5) 



We here consider only the weak-field regime 77$!' <C Siq, 
where only a negligible population is transferred to the 
intermediate levels (n ± 1). In this case and under the 
above laser fields, the transition from |10n) to |een) can 
be anticipated, and the effective Rabi-frequency may be 
evaluated in second-order perturbation theory 
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Em — EiOn 



(6) 



-^int' ^^^'^ intermediate states 

Ee + 

Sio)) and |le(n±l)) 
E, + {n± l)!^ and 



where Hmt 

\m) are |eO(n ± 1)) (with the corresponding E^ 
Eq + {n± l)v and ujm ~ ^ei ± (j^ 
(with the corresponding Em = E ^ -\- tije, 4 
± (J^ - <5io)) (see Fig.3) lll[3. 
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FIG. 3: Schematic diagram of the resonant transition from 
IIOti) \een) for the two ions in a sub-system. The two ions 
are illuminated by two different bichromatic laser fields. Only 
the transitions involving the intermediate states with (n ± 1) 
can occur. 

Using the same scenario, we can obtain the effective 
Rabi-frequencies of the transitions between (|01),|aa)) 
and |ee} respectively as 



(5oi 



(7) 



For simplicity, we may take (Sio = = ^aa = ^ here- 
after. Combining Eqs.(6) and (7) and in the rotat- 
ing frame, the effective Hamiltonian of the sub-system 
(jiij2) with the lasers on can approximately be written 
as 



Heff = [r!i|ee)(10|+17o|ee)(01| 
+na\ee){aa\]+ H.c. 



(8) 
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in the designated operational subspace spanned by {|01), 
10), |ee), |aa)}. At this stage, we consider a smaller 
subspace 



Cj := span{\01), |10}}, 



(9) 



and encode a pair-bit code to construct a logic qubit in 
this subspace: \0)l = |01) and = |10). If we choose 
this logic qubit as computational one. such an encoding 
constitutes the well known DFS scheme that is immune 
from the decoherence induced by the system-environment 
interaction in the form of Z (g) B, where Z = cr^i + ct^^ 
and B is an random bath operator, simply because for 
any state 1-0) G Cj, we have Z\ip) = 0. Also note that, 
we can define the corresponding Pauli operators Rx, Ry, 
Rz in the subspace Cj , which can be expressed as 



Rx 

Ry 

R. 



1 



nX^ X " X ~ "y "y 



(10) 



where the ct^, cr^, cr^(j = Ji, J2) are the Pauli matrices of 
the jth qubit. 

We now address in detail how to realize the holonomic 
computation with a DFS scenario. As is well known, to 
achieve a universal set of quantum gates, we need to con- 
struct two noncommutable single qubit operations and a 
nontrivial two-qubit gate. Here we choose the two single 
logic-qubit gates as Ui = e^'^^^", U2 = e*"^^^", and the 
control phase logic-gate as lA^ = e*'^^'^^''^ similar to 
those discussed in Refs.flllllSj. 

Denoting also \E)l = |ee), \A)l = \aa), the effective 
Hamiltonian of Eq.(8) may be rewritten as 



H,ff = [fli\E)L{l\+n„\E)L{ 
+na\E)L{A\]+H.c.. 



(11) 



To implement the operation Ui, we choose flo — 0, fli — 
51 sin 6*, and fla = flcosOe'"^, where fl is the absolute 
magnitude of the effective Rabi-frequency for the gate 
operations. 9 and ip are the control parameters in the 
parameter space M. In this case the dark states of the 
Hamiltonian(ll) take the form 



1^0) 
1^1) 



|0)l, 
COS 



-lip 



(12) 



We then let ili and ^la change adiabatically from the 
point O {9 = 0) along a close path C in the parameter 
space M, i.e., 9 and ip take a cyclic evolution. Thus for an 
arbitrary initial state 1^*0) in the computational subs pac e 
Cj , the state evolves as |*o) ^ U{C)\^o) El El 111, 
where 



C/(C) = exp ^ A"" 



(13) 



is the non-Abelian holonomy associated with the path 
C and A'^ = A^dXf^ is the U(2)-valued connection 
expressed as 



A^^ = {D,{X)\^\D,{X)) 



(14) 



where {9 or (p) are the coordinates in the parameter 
space. From Eqs.(12)-(14), we can derive the correspond- 
ing connection as 



^^ = -zisin2 (0)(l-i?,) 



(15) 



^ '2 

In this way the holonomic operation is achieved as 

[/(C) = e-*'^ie*«% (16) 

where 



61 = - ® sin^ (6) dif. 



(17) 



Next we deduce the holonomy U2- We parameterize 
Hq ~ rising cos iy9, fli = r2sin6'siniy9, ^la = ^cos6. Then 
the two dark states of Hamiltonian(ll) are given by 

\Do) = cos9cos(p\0)l + cos9sinip\l)L — sin9\A)L, 
\Di) = cos^|1)l -sin^|0)L. (18) 

Taking into account another adiabatic cyclic evolution of 
the parameters 9 and ip, we have 



A^p = —i cosORy. 



Therefore, we obtain 



(19) 



where 02 is the geometric phase factor determined from 
the integral 



^2 = — J) cos 9dip. 



(20) 



So far, we have illustrated the implementation of two 
single-qubit rotations around the y and z axes in the DFS. 
The combination of Ui and U2 allows us to perform any 
single qubit operation. To accomplish a set of univer- 
sal gates, we need to realize one more two-logic-qubit 
controUcd-phase gate, which is usually more crucial and 
important. In our scheme, to realize the two qubit gate 
C/|''' in the DFS, we propose to couple the jith and kith 
(see Fig.l) using the bichromatic laser field, and set the 
j2th and k2th^s energy levels be decoupled. In detail, 
we use two resonant laser fields ujn ~ Uei + {ly — 611) 
and ujii = ujei — {v — Sn) to illuminate the ions ji 
and fci, where u still denotes the phonon frequency of 
ions. We denote the effective resonant Rabi-frequency 

as rij^'^^ The other two pulses UJaa = ^ea + [v — 5aa)^ 

^aa — — {v — 5aa) are uscd to rcsonatc the transitions 
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between \a) and |e) in both the jith and kith ions, with 
the corresponding Rabi-frequency as r^^Q*^^ . As a result, 
denoting \EE)l = |eO)j|eO)fe and \AA)l = \aO)j\aO)k, 
the effective Hamiltonian of the two sub-systems j and k 
with the lasers on may be written in the coded space as 

i/f)^ = nii\EE)L{n\ + nAA\EE)L{AA\ + H.c, (21) 
where 

"11 — "11 — 7 1 

Oil 

Oaa 

From the above Eq. (21), similar to the case in 
the single qubit operation, we have the two zero- 
eigenvalue eigenstates (dark states) as |00)l and 
cos6'|11)l — sm 9e^'^ \AA) L under the parametrization of 
|f2ii|/|£l^^| = tan 6* and (p ~ — (p^^ , where (p"^^ and 
(fi^^ are the phase factors of ^Iaa and flu. When the 
control parameters (0, ip) take a cyclic evolution from 
6' = 0, the component |11)l will accumulate a Berry 



phase While at the same time the other compu- 
tational bases |00)l, \01)l, \10)l are decoupled from 
the Hamiltonian(21). This process corresponds a two- 
qubit controUcd-phasc shift gate U-^'^ = e^'^^^^^^^ (^l in 
the DFS. 

We have proposed a holonomic quantum computa- 
tion scheme in the DFS with trapped ions. We have 
illustrated how to achieve two noncommutablc single 
logic-qubit gates and a two-logic-qubit controUed-phase 
gate using the non-Abelian holonomy. Note that in our 
scheme, the computational bases are all encoded in the 
subspace immune to the cr^-typc of decohcrencc. It is 
believed that the adiabatic holonomies combined with 
the DFS will provide additional virtues against decohcr- 
encc and stochastic errors in controlling operations. It 
seems feasible to realize our scheme in physical systems 
like trapped ions, though it is very challenging experi- 
mentally. 
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